A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
13.07.1, 09.10.1, 08.13.2, 02.19.1 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



Gamma-ray burst afterglows from jetted shocks in wind 
environments 



£2 , L.J. Gou 1 , Z.G. Dai 1 2 , Y.F. Huang 13 , and T. Lu 1 ' 2 

. Department of Astronomy, Nanjing University, Nanjing 210093, P. R. China 

O ■ 2 IHEP, Chinese Academy of Sciences, Beijing 100039, P.R. China 

(N ' 3 



O 



> ■ 

-xf : 

(N ■ 

o ■ 

o : 
o ■ 

Or 

6 ■ 
a: 



Astronomical and Astrophysical Center of East China, Nanjing University, Nanjing 210093, P. R. China 



Received 6 June 2000; accepted 22 September 2000 

Abstract. Gamma-ray bursts with long durations are 
widely thought to arise from the collapse of massive stars, 
where the wind environment is unavoidable. It is also be- 
lieved that 7-ray bursts come from jets. Considering these 
two points in this paper, we calculate the evolution of a 
highly collimated jet that expands in a stellar wind envi- 
ronment and the expected afterglow from such a jet. We 
use a set of refined dynamical equations and a realistic 
lateral speed of the jet, and find: (1) There is no observ- 
able break at the time when the Lorentz factor of the jet 
is equal to the inverse of its initial half-opening angle. (2) 
No obvious break appears at the time when the blast wave 
transits from the relativistic to the non-relativistic phase. 
(3) For the wind case, there is no flattening tendency even 
up to 10 9 s. (4) Compared with the homogeneous medium 
case, our calculated flux is weaker in the stellar wind case. 
Finally, we find that two kinds of GRB models (neutron 
star mergers and massive star collapses) may be discrim- 
inated in our numerical results. 

Key words: Gamma rays: bursts — ISM: jets and out- 
flows — stars: mass loss — shock waves 



1. Introduction 

Since the discovery of the GRB 970228 afterglow by 
BeppoSAX, research on gamma-ray bursts (GRBs) has 
evolved. Now we know that GRBs are one of the most en- 
ergetic phenomena at the cosmological distance. Optical 
afterglows have been observed from about a dozen GRBs 
(Klose 2000), most of which concentrate at the distance 
scale of z *~ 1, corresponding to the luminosity distance 
of about 3.0 Gpc. In addition, some GRBs' host galaxies 
have been discovered. All these discoveries leave no doubt 
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that GRBs, the origins of which had puzzled people since 
their discovery more than 30 years ago, are of cosmological 
origin. 

Whether a jet exists or not in GRBs is a fundamen- 
tal problem. We believe that jets should exist based on 
the following facts: (1) The isotropic energy release per 
GRB is generally in the range of 10 51 — 10 52 ergs. It can 
be explained well by a stellar-mass progenitor. However 
for two GRBs, GRB 990123 (Kulkarni et al. 1999a) and 
GRB 990510 (Harrison et al. 1999), the isotropic energy 
is so enormous that it is difficult to explain it by any stel- 
lar progenitor model, which forces some theorists to de- 
duce that the radiation must be highly collimated in these 
cases. (2) The steepening of some afterglow light curves 
observed at the optical band is argued as evidence that jets 
exist in GRB radiation (Kulkarni et al. 1999a; Harrison et 
al. 1999; Stanek et al. 1999; Huang et al. 2000a). Rhoads 
(1997a, 1999), Sari, Piran & Halpern (1999), Meszaros & 
Rees (1999b) have shown that the lateral expansion of a 
relativistic jet will lead to a more rapid deceleration, caus- 
ing a sharp break in the afterglow light curve. For GRB 
990123, the power law index of the afterglow light curve is 
a = 1.0± 0.03 in the 2 days after the burst. After 2 days, 
there is a sudden steepening in the light curve (Kulkarni et 
al. 1999b). Similarly, for GRB 990510 the power law index 
changes from a = 0.76 to a = 2.4± 0.02 after t = 1.0 day 
(Stanek et al. 1999). Recently, a rapid decay with a = 1.73 
was found in GRB 970228 (Galama et al. 2000). All these 
breaks may be due to the lateral expansion of jets (Huang 
et al. 2000a). (3) The observed radio flare may provide 
an independent and excellent indication of a jet-like ge- 
ometry in GRBs (Harrison et al. 1999; Kulkarni et al. 
1999). As argued by Waxman, Kulkarni and Frail (1998), 
the radio afterglow from a spherical fireball must rise to a 
peak flux on a timescale of a few weeks, but because of lat- 
eral expansion, the radio afterglow from the forward shock 
of a jet must fade a few days after the burst. Therefore, 
the relative faintness of the observed late-time radio cmis- 
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sion implies the existence of a jet (Sari & Piran 1998a). 
GRB 990510 may be a good example: its radio radiation 
began to decline one day after the burst (Harrison et al. 
1999). (4) The polarization observed in the afterglows may 
also be evidence for jets. Gruzinov (1999) has argued that 
optical afterglows from jets can be strongly polarized, in 
principle up to tens of percents, if co- moving magnetic 
fields parallel and perpendicular to the jet have different 
strengths and if we observe the afterglows at a right view- 
ing angle. However, Sari (1999) argued that even if the 
magnetic fields had a well-defined orientation relative to 
the direction of the shock, the polarization is unlikely to 
exceed 20%. Furthermore, taking into account the dynam- 
ics of jets, the polarization first rises to the peak around 
the jet break time and then decays. (5) Observed light 
curves of some GRBs steepen simultaneously at different 
bands. This may be further evidence that there is colli- 
mated ejecta in GRB radiation (Harrison et al. 1999). (6) 
Observational characteristics of some GRBs are similar to 
those of BL Lacs, in which jets are unavoidably involved 
(Paczyhski 1993; Dermer & Chiang 1999; Cheng, Fan & 
Dai 1999). This implies that GRBs may arise from jets. 

Generally, the broken light curves of some afterglows 
are explained to be due to the expansion of jets into a 
homogeneous interstellar medium (Harrison et al. 1999; 
Stanek et al. 1999; Huang et al. 2000a). However, other 
explanations have also been proposed. For example, when 
a spherical fireball evolves in the Wolf-Rayet star wind, the 
light curve can also steepen (Chevalier & Li 1999, 2000; 
Meszaros, Rees & Wijers 1998; Frail et al. 1999a, b; Dai & 
Lu 1998a). Dai & Lu (1999, 2000a, b) suggested that when 
a shock in a dense medium transits from the relativistic 
phase to the non-relativistic phase, a break would occur 
in the light curve. 



Livio & Waxman (1999) recently discussed the evolu- 
tion of a jet in the wind environment and gave an analyt- 
ical result. They argued that at late stages (particularly 
after the break corresponding to 7 = l/9 ), the light curve 
has a flattening tendency. In this paper, we use some re- 
fined equations to describe the evolution of jets. First, in 
the adiabatic case, when blast wave is extremely relativis- 
tic, its dynamical evolution satisfies the Blandford-Mckee 
(1976) solution. But when it reaches the non-relativistic 
phase, it satisfies the Scdov- Taylor self-similar solution. 
However, the conventional dynamical model can not tran- 
sit correctly from the ultra-relativistic phase to the non- 
relativistic phase. This has been stressed by Huang et al. 
(1998a, b, 1999a, b). Here we use the refined dynamical 
equations proposed by Huang et al. (1999a, b, 2000b, c, 
d), which can describe the overall evolution of jets from 
the ultra-relativistic phase to the non-relativistic phase. 
Second, for the lateral expansion speed of jets, it is rea- 
sonable to assume that it is just the co- moving local sound 
speed c s . Usually, one has taken c s = c or c/\/3 (Rhoads 
1997a, b, 1999; Sari, Piran & Halpern 1999), where c is the 
speed of light. In fact we expect c s to vary with time, and 
especially it will by no means be c or c/\/3 when the blast 
wave decelerates into the non-relativistic stage. Huang et 
al. (2000b, c, d) have given the proper lateral expansion 
speed which depends on the bulk speed of the blast wave. 

Based on these considerations, we calculate the evolu- 
tion of jets in the wind environment from the relativistic 
stage to the non-relativistic stage and compare numeri- 
cal results with those of jet evolution in the homogeneous 
medium. We describe our model in section 2. Our detailed 
numerical results are presented in section 3. Section 4 is a 
brief discussion of our final results. 



2. Model 

2.1. Basic Equations 
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dta 
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M ej + em sw + 2(1 - e)jm s 



The most important problem in GRB research is the 
energy mechanism. It is widely believed that GRBs with 
long durations come from the collapse of massive stars 
(Fryer et al. 1998, 1999a, b; In't Zand 1998; Ruffert & 
Janka 1998, 1999; Woosley, MacFadyen & Hcger 1999; 
Rees 1999; Meszaros & Rees 1999a). For a massive star 
at the end of its evolution, it throws away its envelope 
and the core collapses into a compact object, producing a 
jet. Then due to interactions between different shells inside 
the jet and between the jet and its surrounding medium, 
a GRB and its afterglow are produced, respectively. The 
recently discovered connection between supernovae and 
GRBs provides a strong support to such a collapsar model. 
At present, there are three GRBs which are most likely 
connected with supernovae: GRB 980425 (SN 1998bw) 
(Galama et al. 1998; Iwamoto et al. 1998), GRB 970228 dR 
(Galama et al. 2000; Reichart 1999), GRB 980326 (Bloom = @cry (7 + V7 2 - 1J , 

et al. 1999). Recently, two other GRBs were added, i.e. 
GRB 970514 (SN 1997cy) (Germany et al. 2000; Turatto, 
Suzuki & Mazzali 2000) and GRB 980910a (SN 1999e) 
(Kulkarni & Frail 1999c; Thorsett & Hogg 1999). 



The overall evolution of a jet can be described by the 
refined equations (Huang et al. 2000b, c, d) 



x2ttR 2 {1 - cos 6- } )pf3cj{j + ^7 2 - 1), 



= 2ttR 2 (1 - cos 0j) pfley (7 + ^7 2 - l) 



(i) 



(2) 



(3) 



(4) 
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where 7 is the Lorentz factor of the jet, m sw is the total 
mass of the swept-up medium, R is the radius, 0j is the half 
opening angle of the jet, s the observed time, p is the 
mass density of the interstellar medium, M c j is the initial 
mass of the jet, and e is the radiation efficiency of the 
jet. Equations (1) — (4) describe the overall dynamical 
evolution. However before evaluating them numerically, 
we should give the expressions for c s , e, and p. 

2.2. Sound Speed 

The lateral expansion is determined by the co-moving 
sound speed. The simple assumption of c s = c/>/3 is un- 
reasonable in the present paper. Huang et al. (2000b, c, d) 
give the proper sound speed which depends on the bulk 
Lorentz factor. Here we give a brief derivation. Kirk & 
Duffy (1999) have derived 



IP 
P' 



(7 ~ 1)P' 
(7 ~ !)p' + IP 1 



(5) 



where p' and p' are the co-moving mass density and pres- 
sure, respectively, and 7 is the adiabatic index. Dai, Huang 
& Lu (1998, 1999) obtained a simple and useful expres- 
sion for 7, 7 = (47 + 1) / (37). Since e' = 7p'c 2 and 
p' = (7 — l)(e' — p'c 2 ), it is easy to get (Huang et al. 
2000b, c, d) 



1 



W-WT- 1 ) 1 + 7(7-1) 



(6) 



In the ultra- relativistic limit (7 ^> 1,7 w 4/3), equa- 
tion (^) becomes c 2 — c 2 /3; and in the non-relativistic 
limit (7 ~ 1,7 ~ 5/3), we simply get c 2 = 5/3 2 c 2 /9. So, 
equation (^) is a reasonable expression and will be used in 
our model. 



2.3. Radiative Ef, 



iciency 



As usual we assume that the magnetic energy density in 
the co-moving frame is a fraction £3 of the total thermal 
energy density (Dai et al. 1998, 1999) 



B' 2 , 2 77 + 1 , 2 
-^=£B-T3T ( 7- 1)nmpC 



(7) 



Generally, the afterglow comes from synchrotron radia- 
tion of the electrons accelerated behind the shock (Sari 
et al. 1998; Wijers, Rees, Meszaros 1997). The contribu- 
tion of the inverse Compton-scattering emission is always 
neglected, because it is unimportant particularly at late 
times (Waxman 1997). So we only consider synchrotron 
emission from the electrons. Assume that the accelerated 
electrons carry a fraction £ of the proton energy. This 
implies that the minimum Lorentz factor of the random 
motion of electrons in the co-moving frame is 



7e,! 



£e(7 " 1) 



P 



p-l 



1, 



(8) 



where m c is the electron mass, p is the index character- 
izing the power law energy distribution of electrons, and 
m p is the proton mass. Considering only the synchrotron 
radiation, Dai et al. (1998, 1999) derived the radiative effi- 
ciency of the jet 



t'- 1 

"syn 



( 'syn \ bey 
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where t' n 



67rm e c/((7T-B 7c,min) is the synchrotron 
cooling time with or being the Thompson scattering cross 
section, and t' cx — R/^/c is the co- moving frame expan- 
sion time. For the highly radiative expansion, £ c as 1 and 
t' syn <C t' eK we have e w 1. The early evolution of the 
jet might be in this regime. For the adiabatic expansion, 
£ c <C 1 or t' syn ^> t' cx , we get e ~ 0. The late evolution of 
the jet is believed to be in this regime. In this paper, we 
consider the latter case, i.e. adiabatic expansion. 



2.4. Mass Density 

Huang et al. (2000b, c, d) have considered the case that 
a jet evolves in the homogeneous interstellar medium 
(ISM). In this paper, we consider the case that a jet 
expands in the preburst stellar wind. For massive stars, 
particularly Wolf-Rayet stars, the typical wind-loss rate 
is M w 10~ 5 — lO^Afoyr -1 , and their typical speed 
is V w w 1000-25000km ■ s^ 1 (Willis 1991). According to 
4:irR 2 pV w = M, we can easily get 



p = 5.02 x 10~ 18 g • cm" 3 
/ M 



R 



lO^Afoyr- 1 J \ 1000 km -s 



10 15 cm 



(10) 



We will use this kind of mass density structure in the 
paper. 

2.5. Electron Energy Distribution 

In the absence of radiation loss, the distribution of the 
shock accelerated electrons behind the blast wave is usu- 
ally assumed to be a power law function of electron energy, 



dN' 

— — - OC 7 e P , (7o,min < 7c < 7c,max), 
«7c 



(11) 



where 7 G ,max is the maximum Lorentz factor, 7 e ,max = 
lO^B'/lG)- 1 / 2 (Dai et al. 1998, 1999), and p usually 
varies between 2 and 3. However, radiation loss may play 
an important role in the process. Electrons with differ- 
ent Lorentz factors have different radiation efficiencies. 
Sari, Piran & Narayan (1998) have derived an equation 
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Fig. 1. Evolution of the Lorentz factor. We have taken the 
"standard" parameters: E /Q, = 1 x 10 54 ergs/47r, 70 = 300 



(i.e. M c j/fto 



O.OO2M /4tt), 



0.02, p 



2.5, 



£ e = 0.1,6*0 = 0.2. The solid line corresponds to the wind 
environment (M = 1O~ 5 M yr -1 and V w = 10 s kms -1 ), and 
the dashed line to the homogeneous ISM (n = 1cm -3 ). The 
time that 7 = 2 is about 6 x 10 6 s in the ISM case, and 10 8 s 
in the wind case. 



for the critical electron Lorentz factor, j c , above which 
synchrotron radiation is significant, 



7c 



6nm c c 



(12) 



Electrons with Lorentz factors below j c are adiabatic, and 
electrons above 7 C are highly radiative. 

In the presence of a steady injection of electrons accel- 
erated by the shock, the distribution of radiative electrons 
becomes another power law function with an index of p+ 1 
(Rybicki & Lightman 1979), but the distribution of adia- 
batic electrons is unchanged. Then the actual distribution 
should be given according to the following cases (Dai et 
al. 1998d, 1999): 

(1) For 7 C < 7 C 

.min i 

dNL 



Ci = 



= Cl 7e 



(p+1) 



7e 



/c,max 



(7e,min ^ 7c — 7c, max) 5 (-^) 



N clc , (14) 



where N e \ e is the number of radiating electrons in a 
ring between 9 and 9 + dd with 9 being the angle be- 
tween the velocity of emitting material and the line of 
sight. 

(2) For 7c,min < 7c < 7e,max, 

dK _ ( C 2lc P , (7c,min < 7c < 7c), 

d Jc lC 3 7c" (P+1) , (7c < 7c < 7c,max), 

where 

C 2 = C 3 / 7 c , 



(15) 



'c,min lc ^ ic /c.max 



7c (p- 1) 
(3) For 7 c > 7c 



P 



N e l e . 



dK 

djc 



CAl- p , (7c ,min < 7c < 7c,max), 



(17) 



(18) 




o 

log t(s) 

Fig. 2. Evolution of the shock radius (R). Parameters and line 
styles are the same as in Fig.l. 



where 



Ci = — 



p-1 



■y~ p -■y 1 ~ p 

'c.min /c.max 



N ele . 



(19) 



2. 6. Formulae of Synchrotron Spectrum 



In the co-moving frame, the synchrotron radiation power 
at frequency v' from electrons is given by (Rybicki & 
Lightman 1979) 



m e c^ J 7o min V 07e / V !-. 



where e is electron charge, v' c = 3j^eB' /(47rm c c), and 



r-too 

F(x) = x K 5/3 (k)dk, 

J X 



(21) 



with K 5 / 3 (k) being the Bessel function. We assume that 
this power is radiated isotropically in the comoving frame, 



dP'{y') P'{v') 



dfi' 



47T 



(22) 



Defining /i = cos 9, we can derive the differential power 
in the observer's frame (Rybicki & Lightman 1979; Huang 
et al. 2000b, d), 



dP(v) 



1 



dP'{v') 



1 



P'{y') 



(16) dQ 7 3 (l-/3^) 3 rfO' 7 3 (1-/3 M ) 3 4tt 



,(23) 
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7(1 - nP) ■ 



(24) 



Here quantities with prime are measured in the comoving 
frame, and quantities without prime are in the observer's 
frame. Then the observed flux density at frequency v at 
certain angle is 




4 

log t(s) 

Fig. 3. Evolution of the half opening angle(#). Parameters and 
line styles are the same as in Fig.l. 



S v {n) = -r- 



dPjy) A s 
dfl Dl 



7 3 (1 - [3^if 4ttL>2 



P' ( 7 (1 - tf)v) , 



(25) 



where A s is the area of the detector and Dl is the lu- 
minosity distance. The observed flux density at a given 
frequency is obtained by integrating over the shock front 
within the jet boundary 6 y 

3. Numerical Results 

In our model, we use the following initial values or pa- 
rameters as a set of "standard" parameters: Eq/VIq = 
1 x 10 54 ergs/47r,7o = 300 (i.e. M cj /ft = O.OO2M /4tt), 
£| = 0.02, p = 2.5, Co = O.l,0o = 0.2. For simplicity, 
we assume that the expansion during the whole stage is 
adiabatic, i.e. e = 0. 

Figure 1 shows the evolution of the Lorentz factor. We 
see that the bulk Lorentz factor changes very slowly at late 
times in the wind case, compared with that in the homo- 
geneous ISM case. For the wind case, using the standard 
parameters, we see that at 10 8 s, the blast wave deceler- 
ates into the non-relativistic stage (here we let 7 = 2 be 
the critical point between relativistic and non-relativistic 
phases). In the homogeneous ISM case, the blast wave 
evolves into the non-relativistic stage at about 10 7 s. Fig- 
ure 2 illustrates the time dependence of the shock radius. 




Fig. 4. R band afterglow light curves. The thick solid line 
corresponds to a jet with "standard" parameters, and viewing 
angle 6> obs =0. 



In figure 3 we present the evolution of the jet open- 
ing angle. During the ultra-relativistic phase the angle in- 
creases only slightly, because at this time, lateral expan- 
sion speed can be thoroughly neglected when compared 
with the blast wave speed itself. However at the Newto- 
nian stage, 9 increases quickly from 0.6 to 1.4. 

In figure 4, we show the R-band light curve. Please note 
that at the end point of each curve the average electron 
Lorentz factor is already as small as 7 , m in = 5, corre- 
sponding to a bulk Lorentz factor of 7 = 1.05, which im- 
plies that the jet is completely in the Newtonian regime. 
Rhoads (1999) as well as Livio & Waxman (1999) has pre- 
dicted that the light curve will show a break when the bulk 
Lorentz factor is 7 w 1/0 O (here the blast wave is still at 
the relativistic stage ). In figure 5, we give the evolution 
of the time index of the afterglow. We expect that the 
break in the wind case is less obvious than that found by 
Huang et al. (2000a, b, c, d) in the homogeneous ISM. 
Our numerical results verify this expectation. 

Combining figure 4 with figure 1, we can see that 
both in the wind environment and in the homogeneous 
ISM environment there is no observable break during the 
relativistic stage, which is consistent with the results of 
Panaitescu & Meszaros (1998), Moderski et al. (2000), 
Huang et al. (2000a, b, c, d), and Wei & Lu (2000). We 
expected that there would be an obvious break during the 
trans-relativistic stage, i.e. transition from the relativistic 
stage to the non-relativistic stage. We don't find such a 
break but rather a smooth curve, which can be seen clearly 
in figure 5. In a uniform density medium the increase of 
the index in the power-law of the light curve is 1.07 during 
about two and a half decades in time. For a pre-ejected 
stellar wind a increases by 0.8 over 5 decades. Therefore, 
as argued by Kumar & Panaitescu (2000), a break in the 
light curve for a jet in a wind model is unlikely to be 
detected. 
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4. Discussion 

Following Huang et al. (2000a, b, c, d), who considered 
the evolution of a jet in a homogeneous ISM environment, 
we investigated the detailed dynamical evolution of jets 
and their afterglows for the wind case from the ultra- 
relativistic stage to the non-relativistic stage. Recently, 
Kumar & Panaitescu (2000) considered the evolution of a 
jet in stratified media. Compared with their studies, our 
model is refined in the following aspects: (1) Kumar & 
Panaitescu (2000) considered the evolution of an adiabatic 
jet, while we study the expansion of a partially radiative 
realistic jet. Furthermore, the dynamics presented are ap- 
plicable to both ultra-relativistic and Newtonian jets, so 
we could follow the overall evolution of a jet using a set of 
differential equations. (2) Similarly, Kumar & Panaitescu 
(2000) did not consider the variation of sound velocity 
with time. We describe the lateral expansion of jets with 
a refined and more reasonable sound speed expression, 
which varies with the bulk Lorentz factor (Huang et al. 
2000b, d). 
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-2.1 



-i 1 r - 




-0.9 



■1.2 



■1.5 



-2.1 
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Fig. 5. The evolution of the time index of the afterglow light 
curve at the R band. Parameters and line styles are the same 
as in Fig.l. 



In addition, we considered the evolution of electron 
distribution with time. Despite these differences, one of 
our results is similar to that of Kumar & Panaitescu 



2. The transition from the ultra-relativistic phase to 
the non-relativistic phase is also very smooth; the ex- 
pected obvious break does not appear. This is very simi- 
lar to the behaviour of isotropic fireballs (Wijers, Rees, & 
Meszarosl997, Huang et al. 2000b). 

3. If we use the same parameters (except for the dif- 
ferences in number density for the wind case and for the 
homogeneous ISM), we find that the flux density in the 
wind case is obviously weaker than that in the homoge- 
neous ISM case. This property is consistent with Chevalier 
& Li's (2000) result. 

Two currently popular models for GRB progenitors 
are the mergers of compact objects (neutron stars or black 
holes) and the explosions of massive stars. It is widely be- 
lieved that GRBs produced by the former model occur in 
the ISM with density n ~ lcm -3 and GRBs produced by 
the latter model occur in the preburst stellar wind enviro- 
ment with mass density p cx R~ 2 . As argued by Chevalier 
k Li (2000) and Livio & Waxman (1999), both ISM and 
wind cases should show the same emission feature dur- 
ing the lateral spreading phase, and in particular on a 
timescale of days, the wind density is similar to typical 
ISM densities so that an interaction with the wind would 
give results that are not different from the ISM case. If 
GRBs are beamed, thus, their optical afterglow emission 
could not be used to discriminate the massive prognitor 
model from the compact binary progenitor model. How- 
ever, our numerical results show that their optical after- 
glow emissions arc different, particularly several days after 
the burst. Thus, it may be used to discriminate the two 
models from each other, and further observations may ver- 
ify our numerical results. 
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in the light curve. 
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